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Abstract. 

Pure gauge lattice QCD at arbitrary D is considered. Exact integration over link 
variables in an arbitrary D- volume leads naturally to an appearance of a set of surfaces 
filling the volume and gives an exact expression for functional of their boundaries. The 
interaction between each two surfaces is proportional to their common area and is real- 
ized by a non-local matrix differential operator acting on their boundaries. The surface 
self-interaction is given by the QCD 2 functional of boundary. Partition functions and 
observables (Wilson loop averages) are written as an averages over all configurations 
of an integer-valued field living on a surfaces. 



1 Introduction. 

Pure gauge QCD (quantum gluodynamics) in D dimensions has been reformulated by 
K.Wilson in terms of collective (lattice) variables [0: 

S=y £tr(tf / + t# (1) 
A G f 

where sum goes over all faces / of D-dimensional lattice, Uf — ILe/ U{ with [/, being 
unitary matrix (U(N) or SU(N)) attached to i-th link of a lattice, A Q is the bar (lattice) 
coupling constant. In continuum limit, at any D < 4, A goes to zero. At D = 2 
A ~ e 2 , at D = 3 A D ~ e and at D = 4 A Q ~ — where e is the linear size of a 
lattice. 

The partition function is defined as an integral over all link variables, 



l[dU t e s , (2) 

i 

while the simplest observable, the Wilson loop average, is 

W(C) = -^- fUd^ehrdiU,) , (3) 

iVZ/ J i jec 

where C is a closed contour. 

In order to integrate exactly over unitary matrices one makes the group Fourier 
transformation of the face variables, first manifestly used in D = 2 by A.A.Migdal [Q: 

e ^f^ =^d r A r (^) Xr (U f ) , (4) 



d r A r (£) = / dU Xr (U) e ^+^ , (5) 

where r is irreducible representation of the gauge group, Xr{U) is its character and 
d r = XrOO is its dimension, 

dr= II (1 + ^^) • (6) 



l<i<j<AT 



Here and below we use the standard parametrization of r by its highest weight com- 
ponents, ni > ... > n N , associated with a lengths of lines in the Young table. 
Direct calculation of (H) gives 



. /AN 1 , T f 2N^ 

r \Y) = T det « J "i-*+nirJ ' ^ ^ 



■A ^ d r lJIL ^ J ^\ 
where I n (x) is the modified Bessel function. 
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In D = 2, an exact solution of the model can be obtained using only orthogonality 
of characters, 



fdU X r 1 (AU)xr a (rtB) = 6 ri 



Xn(AB) 



''2 



d 



2> 



The result of integration inside a disk of area A (therefore, its lattice area is A/e 
gives functional of boundary [] T M: 

Z^ t {T)=Y.dX le2 Xr{T). (9) 

r 

In the continuum limit, e — > and A D = Ae 2 , we need an asymptotic expansion of A r , 
which can be found from (H) by the saddle point method to give: 

a, - i - ^f 1 + o(M) (10) 

up to representation- independent factor. Here, C 2 (r) is the eigenvalue of the second 
Casimir operator, 

N 

C 2 (r)=Y / n l (n t + N + l-2i) . (11) 

8=1 

Thus, the continuum limit of (Q) is defined by the substitution 

- exp ( - ^C 2 (r) 
and takes the form 

Z(T) = 2>exp ( - A-C 2 ( r )) Xr (T) . (12) 

The same can be done in the case of non-trivial topology. The results are || Q 

(i) functional of boundaries Tj, % = l,...,n of a sphere with n holes and of area A 
(continuum coupling constant A is absorbed into the area): 

A n 

Z(T U ...,T n ) = Y,d 2 r - n exp ( - ^C 2 (r)) nx,(r,) , (13) 

(ii) partition function for a closed surface of a genus g and of area A: 

Z a (A) = Y. d l (1 - 9) ^p{-4jC 2 (r)) , (14) 

(iii) Wilson loop average: 

m a 

w 9 (c)= y: *r,, m ii< Hi)ex p(-^w) > (is) 

n,...,r m k=l /iV 



*By "boundary", here and below, we equally imply the geometrical boundary and the product of 
unitary matrices attached to it. 
2 See also 01. 
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where m is the number of windows, A k is the area of a window, g k is the "genus per 
window" and coefficient $ n ...r m is the U(N) (SU(N)) group factor dependent on the 
contour topology (see Ref.0 for details). 

In D > 2, the orthogonality condition (§) is not enough to perform an integration 
over link variables since there are more than two plaquettes match at each link. For- 
mally, we still could integrate using known formulas for tensor product of irreducible 
representations. This results in a sums over internal spaces of representations. The 
Clebsch-Gordan coefficients entering in these expansions are not known in any com- 
pact and general form. Besides, after integration we should perform a summation of 
a resulting expressions weighted with these Clebsch-Gordan coefficients, which makes 
the problem to be extremely difficult. However, the problem becomes less hopeless if 
one guesses that the three steps, namely, expansion into representations, integration 
over link variables and then back re-summation, could be performed in one step if some 
adequate variables are found, which should be a combination of one-link integral and 
of a sum over representations. 

In the next two sections we present such a variables and describe the procedure of 
integration in an arbitrary (lattice) D-volume. 



2 One-link integral. 

Though our final results will not depend on choice of a lattice, we start for concreteness 
from the quadrilated regular lattice. 

There are 2D — 2 plaquettes interacting through each link on a quadrilated D- 
dimensional lattice. The one-link integral has the form 

2D-2 

JV 



dU n e£ trW ^ +h - c - , (16) 

k=l 



where by W k we denote a product of three other unitary matrices in a k-th plaquette 
(see left side of Fig.p] for 3D example). In the heat kernel framework, fll6|) is equivalent 
to 



2D-2 , 2D-2 



E II < e-&°>™ fdUjl Xn(W k U) . (17) 

f^T^ ^ 7 1 7 1 



n...T2D-2 k=l k=l 



Expressions (|T6| ) and ( |iT|) are identical up to the 0(X ) order. The difference appears 
in order 0(\ 2 ). 

Making shift U WjlJ (the Haar measure dU is invariant with respect to such 
transformations) we see that the integral fll6|), QlTD depends on 2D — 1 boundaries 
W k W], j + k. The shift can be made in 2D — 2 possible directions corresponding 
to 2D — 2 possible choices of j. The picture corresponding to one of a four possible 
directions of gluing in D = 3 is presented on Fig.|l|. 
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Figure 1: Gluing in three-dimensions. 



Thus, after one-link integration we obtain 2D — 1 two-plaquette surfaces with some 
interaction between them. We are going now to calculate this interaction and to con- 
tinue the procedure consequently for all links. 

Let us start from the heat kernel representation (|TT| ) and consider for simplicity the 
case of only three plaquettes, 

£ d ri d r2 d r3 e -^(n)+c 2 (r 2 )+c 2 (H0) f dUx ri (WiU)xr 2 (W 2 U)xr 3 (U) , (18) 

where W\ and W 2 are the boundaries of the two-plaquette surfaces. 
We are going now to calculate exactly the quantity 

Y,d r e~^ c ^ I dU Xri (WlU) X r 2 (W2U)Xr(U) . (19) 
r 

A direct strategy would be to integrate first over U and then to take the sum over 
r and over its internal sub-space. As we have mentioned in Introduction this way is 
extremely difficult. Instead of that, we first replace the sum over r by the original 
Wilson term exp ^-tr (U + W) and then derive the heat kernel exponent as a first 
order non-zero term in A G -expansion of an integral over U. This will give an adequate 
variables for (|T9"D. 

To derive A Q -expansion of the integral 

dU ^V+v\ ri{Wl U) X r>{W 2 U) , /(Ao)= fdUefrV^ , (20) 



/(A ) 



we diagonalize U, U = QuQ^ where u = diag{ e*^ 1 ,..., e l< ^ N }, and integrate over 
diagonal u near the saddle point u = I. Then, ( |2"0"| ) takes the form: 

Xn{Wi)Xr.{W 2 )(l - ^[C 2 (n) + C 2 (r 2 ) - 2C 1 (r 1 )C 1 (r 2 )]) 

(21) 

+ X o— ±- E / dto-jrr- XrAWiVutf)— xr.(w 2 nurt) + o(\ 2 ) , 

N fc^J 0(p k u=I d<p k u=I 
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where the relation 



U=I 



Xr(WttuQ)) = iCi(r)x r {W) 



iV 



Ci( r ) = E Uk 

k=l 



(22) 



has been used (Ci(r) is the first Casimir operator eigenvalue). Indeed, 

JIM > 



N p) N p, 

E arl=M wnun ^ = Y.TL(rtwn) £ ( 



U=7 



(23) 



k=l "V* " ' a k=l 

where rig is the matrix element of an r-th irreducible representation. There are only 
diagonal matrix elements in (^) since u is diagonal. The matrix label a is parametrized 
by the following Gelfand-Zetlin patterns H: 



a 



rv 



\ 



n 



n 



JV-l 
N-l 



(24) 



where 



ni > ni 1 > ni 



<>k — l0 k — "fc+i > 

(the top level numbers are coincide with the highest weight components: = rik, k 
1,...,N). Then 



(25) 



N 

n 



where c,- is a sum of numbers of j-th level, 



and 5c j = Cj — Cj_i. Besides, 



fe=i 



TV 

E = c " 

k=l 



c = 



Ci(r) • 



(26) 



(27) 



(28) 



Then, the relation (|22| ) becomes obvious. 

To make further calculation we use the formula ( Weyl second formula for charac- 
ters) : 

Xr {W) = det i3 a nt . i+j (W) , (29) 

where numbers rij's are the highest weight components of r and a n (W) is the character 
of representation {n, 0, 0} Q This character can be written as 



a n (W) = w n>1 (tr W) n + E «>n, g (tr WT^tr W q 

q=2 



(30) 



3 In formula (p9|), <7„ with n < might occur. The convention for that case is a n = 0. 
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(an exact form of a coefficients w n ^ is not needed for the following consideration). Since 

a n 

a n (WQu^) = (Q^WQ) kk a n ^(W) + ^ qw q (Q^W q Q) kk (tvW) n ~ 9 (31) 

g=2 



— I- 



and 



we find: 



N 



£ I rfO(O t ^iQ) fcifc (Q t W 2 n) fefe = jJ-^(trW 1 trW 2 + trW 1 W 2 



k=l 



(32) 



(N + 1) £ / dfi^r ^(^JM*)— - <7 ro (W 2 fiutf) 



(33) 



^tr 9^ tr d W2 + tr cV x cV 2 )a n (W 1 )a m (W 2 ) 
where matrix elements of dw are defined by 



TV 











dW k 



i.e., matrix elements Wj k of fundamental representation has to be considered as an 
independent variables (in other words, an action of derivative (|34]) is defined on the 
GL(N) group). An important property of derivative dw is its invariance under right 
group transformations [], 

d w = dwv , V e GL(N) . (35) 



-Wj m — 5j k 5 r , 



(34) 



It is not difficult now to generalize (|33|) to the case of arbitrary representations: 
= (tr 9^ tr <9jy 2 + tr d Wl d W2 ) Xn ( Wi) Xr 2 ( W 2 ) • 

Finally, since 

tvd wX r(W) = iC 1 (r) X r(W) 
(compare with (^)), equation ( p0[ ) takes the form Q: 

= i 1 - ^[^(ri) + C 2 (r 2 ) - 2tid Wl d W2 ])xr 1 {W 1 )Xr 2 iW 2 ) + 0(\ 2 ) . 



(36) 



(37) 



(38) 



4 We could define dw as 



N 

(V) =«E^ 



Then, 5^ is invariant under the left shift — > FVF. 

5 In what follows we put /(A ) = 1 since this factor is representation-independent and, therefore, 
precisely cancels in any physical quantity (normalized to partition function). 
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Therefore, 

Y,d r e-^ c ^ J dU Xn (W 1 U)x n (W 2 U)xr{U) 



(39) 



e ~^(C 2 (r 1 )+C 2 (r 2 )-2tre Wl ^ 2 )^^ H/i ^ r ^ w/ ^ _ 



As explained above, the result (|39|) is exact, i.e., valid for arbitrary A Q . It can be 
also viewed as follows. Equation ( Jl9| ) can be written as 

£d r e^ A ^ / dU Xri (WiU)Xr 2 (W2U)Xr(WU) 

r J 

(40) 

= p T^Nf^ 1 ! dU ^{WxU) Xr2 {W 2 U)Y,d rX r{WU) , 

where Aw is an invariant Laplace-Beltrami operator which is the differential operator 
with respect to parameters of W. A character of an arbitrary irreducible representation 
is the eigenfunction of A, 

A wXr {W) = -C 2 {r) Xr (W) . (41) 
Using the completeness condition for a characters, 

J2d rXr (W)=S(W,I) , (42) 



we write ([|0|) as 

oo \ k 

£ 777^ A w=i Xn {W 1 W) Xr2 (W 2 W) . (43) 



Actually, we have checked above that 



&w Xri (WiW) Xr2 (W 2 W) = -(C 2 ( ri ) + C 2 (r 2 ) - 2trd Wl d W2 ) Xri (W 1 W) Xr2 (W 2 W) . 

(44) 

Therefore, 

^wXr 1 (W 1 W) X r a (W 2 W) = n\C 2 (r 1 ) + C 2 (r 2 )-2trd Wl d W2 ) k Xri (W 1 W) Xr2 (W 2 W) , 

(45) 

which confirms fl3"9"|). 

A generalization to the case of arbitrary number P of plaquettes joining at one link 
is straightforward: 

J2d r e-w^M / dU Xr (U) n Xr P (W p U) 

r J p=l 

(46) 

= ex P-^(E^W-2 E ^d Wp d Wq )l[xr p (W p ). 

ZiV p=l <pcj> p=l 
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Operator exp ^tr dw i dw j acts on a characters in a simple way. In Appendix we 
give several examples corresponding to several low- dimensional representations. It is 
worth mentioning, that although the formula (|29|) is written for non-negative signatures 
(n>N > 0) only, the results ([46|) are valid for an arbitrary representations. 



3 Functional of boundaries and set E. 

Using property ([35|) we proceed further link by link and integrate out all link variables 
inside an arbitrary D-volume. In this way we obtain some set E of (interacting) 
surfaces <Sj. Suppose, all of them are disks (it is clear, at any D it can be done, at 
least in a small enough volume). Let Tj is the boundary of i-th disk and Ai = A(Si) 
is its (lattice) area. To each «Sj corresponds the sum over irreducible representations 
rj. Corresponding functional of boundaries Zj^(^{T}^j has the form of a sum over all 
configurations {r}: 

^({r}) = E e" & «« II d nX ri (Ti) , (47) 

{r} iGS 

^(W, W) - ^ E ^trft^ , (48) 

ieE <ij> 

where Aij = A(SiC\Sj). Formula (|47|) generalizes the expression (12) for the functional 
of boundary in D = 2. The latter corresponds to ( WJf) with E containing only one 
surface. 

There is an infinite number of equivalent sets E. Actually, each set is defined by the 
local gauge fixing, i.e., all such a sets are gauge equivalent. It is clear, observables are 
independent on the choice which should be dictated just by convenience of calculations. 
In D = 3, for example, it is possible to chose E containing surfaces of only disk topology. 
The example is drawn (in projection) in Fig. 2 (a). In this case, all disks are compressed 
and has the form of a closed from one side cylinders of one plaquette width, which fill 
densely a 3D volume. Choosing another direction of gluing at (at least) one link (see 
Fig.l) we obtain another set £' which differ from E not only by smooth deformation 
of a surfaces but also by an appearance of a compact surfaces (see Fig. 2(b)). 

To see what happens when the compact surfaces appear in the set we, first, consider 
the situation when the surfaces with a boundaries are not only disks (1-holed spheres) 
but also spheres with an arbitrary number of holes. It is not difficult to see that in 
such a case we have to replace an expression d ri x n (^i) which appears in ( f4"TD and 
corresponds to z-th disk by d 2 ~ ni llfeLi XrS^i,k) where is the boundary of k-th hole 
{k = 1, ...,Ui) on 2-th surface. 

Thus, the most general expression for functional of boundaries is: 

m 

^({r}) = E e-*«'M*» II II XnOV) • (49) 

{r} i6S k=l 



S 



(a) 



(b) 



Figure 2: Two examples of gauge equivalent sets of surfaces (projection from D = 3). This equally 
can be considered as a set of curves for the Principal Chiral Field model. 

A compactification of a surfaces will become clear in the next section, where we will 
consider a partition functions. 

4 Evolution operator, partition functions and loop 
averages. 

The expression fl49|) can be equally written as 

Z^{{T})=R\{Z{T t ) ) (50) 

R = exp ^ A ij tid Vi d Vj , (51) 

<ij> 

where the 2D boundary functional Z(Ti) is given by ( fi~2|) in the case of a disk or by (|13|) 
in a general case (then, Tj means the set of boundaries {T^i, corresponding 
to i-th surface). Thus, the differential operator J? describes an evolution of QCD 2 in 
QCD D . 

The expression for functional of boundaries ([49"D can serve as the building block in 
construction of partition functions and observables on arbitrary D-manifolds including 
those of non-trivial topology. 

It is clear, the partition function can be obtained from fl49|) by putting all boundaries 
Ti be equal to /, i.e., 

m 

= £ e-^},{3r = ,}) jj d 2- ni jj ^ {v . k) (52) 

{r} iSS fe=l 

where ^({r}, {<9r}) is defined in ([48]). This is a general expression for the partition 
function of U(N) (and SU(N)) lattice quantum gauge theory in D-dimensions. 

To write the Wilson loop average ([3D in the same terms let us consider an arbitrary 
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surface Sc such that C = dSc- Then, the result is 



1 i U i 

W(C) = -f-J2 e-*«-M*-»n e^ Scn5 ^ tra -^-trCn^r ni nXn(r, fc ) - 

S {r} 3 fe=l 

(53) 

Differentiation with respect to C is easy to perform and we obtain an equivalent 
expression: 

W(C) = -^E e-*»««-W*-»tr(n e^V nCIlA) • (54) 
s H j i£E fc=i 

Thus, the loop average takes the form of an average of differential operator 

ltr(n e ^ AiScnS ^) . (55) 

3 

It is not difficult to check that the result is independent on a choice of Sc- 



5 Discussion. 

We represented pure gauge lattice QCD^ as a statistical model of integer- valued scalar 
fields living on a set of interacting surfaces. Apparently, modulo some possible (and 
hopefully fruitful) re-writing, the expressions derived here for the functional of bound- 
aries P^D, for the partition function (^) and for loop averages (|53"D cannot be simplified 
further, unless continuum limit is taken. This can be seen already in the U(l) case 
(QED), where partition function takes the simple form: 

= £ e-*>«»» , 
{«} 

A (56) 
Ss( W) = yE An 2 i - A J2 A ij n i n j ■ 
«es <ij> 

By "simplification" we mean an essential reduction of configuration space {r}. Such 
a simplification is expected in continuum limit. We hope that the derived integrated 
version of the model is a better starting point for analytical study the continuum limit 
than the original, non-integrated version. The problem might perhaps be solvable by 
methods of elementary combinatorics, where the only subtlety is the calculation of the 
surface entropy factor which manifests itself in a change of the set S under refinement 
of the lattice Q. 

The method of integration over unitary matrices and the formula (|3^) (or, in gen- 
eral, (f46|)) can be applied without any changes to the Principal Chiral Field (PCF) 

6 I wish to thank M.Karlincr for discussion of this point. 
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model (for recent developments in the model and for references see 0). The set S 
in this case becomes a set of curves (see Fig. 2) and the partition function is again of 



the form (52), after we replace areas of a surfaces by lengths of curves. Apparently, 
due to its respective simplicity, the PCF model might become the first case where the 
continuum limit can be taken. 

The problem of integration over unitary matrices in QCD is similar to one in matrix 
models of 2d quantum gravity embedded in D > 1 (the simplest case where the problem 
appears is the closed D = 1 target space, i.e., closed chain of hermitian matrices). The 
present method can be applied, after some modification, to these models as well. 

Among possible direct continuations of our analysis let us mention the following: 

(i) Expression fl4"9|) looks suitable for 1/A-expansion. The latter has been recently 
elaborated in a detail for D = 2 @-@ []. It is tempting to apply this technique to 
D > 2. A straightforward strategy would be to find out a "stringy" interpretation 
for the evolution operator flSip . Then the whole theory could be considered as a set 
of stringy models (QCD 2 ), interacting through the operator R. However, a less naive 
strategy is possible, if one makes the 1/ A-expansion of whole expression (^9|), together 
with i?-operator. Then, each term of 1/ A-expansion will take a form of a sum over all 
surfaces from the set E. This requires some new technique, especially with respect to 
i?-operator. 

(ii) The derived exact expression for functional of boundaries (^9|), especially written in 
the form (|50|) , can serve as a starting point for a mean field analysis of the model 0. An 
equation for eigenvalues of i?-operator, under mean field assumption, seems solvable, 
at least on a finite lattices. 

(iii) It is not difficult to recognize that the expression ([52]) for QCD partition function 
can be interpreted at infinite A as the constrained matrix model. The corresponding 
technique was worked out for D = 2 in [|12|]-[|14|]. It would be interesting to solve the 



large N saddle point equation and to see if there is a region of the coupling constant 
where unitary constraint can be ignored, which would indicate an existence of large N 
(apparently, third order) phase transition [13 . 
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Appendix. 



In this appendix, we will demonstrate on a few examples an action of operator 
exp'ytr OaOb (7 is a parameter) on the function Xn(A)Xr 2 (B) ■ Thus, we calculate 

e^ 9 ^ X rM)Xr 2 (B) . (57) 

(a) If ri or r 2 is trivial representation {1, 0, 0}, then (|57]) is equal to 1. 

(b) The first non-trivial example is r\ = r 2 = {1, 0, 0} (fundamental representa- 
tion). In this case we have 

e7 tra A 9 str A tr B = y _L_ tr iq a q bXiI AtiB = tiAtiB cosh 7 - tr AB sinh 7 (58) 

since the even-order derivatives result in ti A tr B and odd-order ones result in trAB. 

(c) n = {1, 0, 0}, r 2 = {1, 1, 0, 0}. Then 

3 27 - e" 7 „ , , „ „ 9N e 27 + 2 e" 7 



P 2 7 — p-7 p 2 7 _|_ 9 p-7 

(59) 

(tiBtrAB - trAB 2 ) + tiA(ti 2 B -tr£ 2 ^ 



3 x 6 
(d) n = {1, 0, 0}, r 2 = {2, 0, 0}. Then 

p-27 _ e 7 p-27 _L 2 e 7 ( 60 ) 

= (trStrAB + trAB 2 ) + tiA(ti 2 B +trB 2 ) . 

3 6 

We leave it to reader to substitute these results in the formula (|39"D to see how it 
works in these cases. 
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